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PREFACE 


Research  under  the  contract  has  been  conducted  to  the  subject 
of  data  compression  for  microwave  imagery.  A  theory  of  Bit 
Compression  has  been  developed  for  compressing  the  aperture  data 
in  a  microwave  antenna  array.  Two  papers  have  been  sent  to  IEEE 
for  publication.  The  first  paper:  "Optimum  Data  Quantization  in 
Microwave  Imagery  and  Its  Effect  Upon  Image  Quality"  is  to  be 
published  in  IEEE  Transactions  on  Acoustics  Speech  Signal 
Processing  of  December  1988,  and  the  second  paper: "An  Exact 
Solution  of  the  Analytic  Equation  of  Image  Quality  from  Optimum 
Quantization  of  Microwave  Imaging  Data"  is  currently  under  review. 

In  this  report,  the  objective  of  this  research  is  reviewed, 
the  major  procedure  of  theoretical  development  is  demonstrated  and 
the  important  results  are  discussed. 

The  mathematical  derivation  will  not  be  shown  in  every  detail 
in  the  report  and  could  be  found  in  the  papers. 
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ABSTRACT 


In  microwave  imaging/  aperture  data  are  digitized/  while 
system  performance  is  determined  by  the  quality  of  the  image, 
which  is  usually  a  Fourier  Transform  of  the  aperture  data.  This 
research  develops  a  methodology  for  finding  a  quantizer  in  the 
aperture  domain  which  optimizes  an  image  plane  quality  criterion. 

For  amplitude-phase  quantization  of  complex  data,  a  novel 
optimality  criterion  is  posed  which  leads  to  a  new  and  simpler 
optimal  quantizer  design. 

The  effectiveness  of  the  theory  is  evaluated  with  computer 
simulation  and  on  experimental  data.  The  theory  explains  why 
fewer  bits  are  satisfactory  in  reconstructed  image,  and  it 
provides  an  analytical  basis  and  a  rule  of  thumb  for  bit 
allocation  between  amplitude  and  phase. 

The  theory  imposes  substantially  no  constraints  upon  the 
distributions  of  amplitude  and  phase  in  the  microwave  data.  It  is 
also  shown  that  the  optimum  quantization  design  derived  from  the 
theory  is  nearly  scene-independent  and  may  achieve  real-time 
performance . 

Finally,  the  theory  provides  useful  guidelines  to  analyze  the 
nature  of  distortion  due  to  data  compression. 


1.  INTRODUCTION 


A  high  resolution  microwave  imaging  system  must  process  huge 
amounts  of  data.  Because  of  the  limitations  on  processing  time, 
memory  space  and  the  cost  of  processor  hardware,  it  is  important 
to  develop  techniques  to  reduce  the  data  handling  requirement. 

The  data  handling  rate  denoted  by  R  is  the  product  of  the 
number  of  the  samples  denoted  by  N  (which  is  the  number  of  antenna 
array  elements),  the  bandwidth  of  the  signal  denoted  by  B  and  the 
bit  per  sample  denoted  by  b  (bit  is  used  as  the  information 
carrying  unit).  Namely, 

R  *  NBb  (bit/second) 

High  cross-range  resolution  requires  large  number  of  elements 
and  high  range  resolution  requires  large  bandwidth.  Therefore, 
for  high  resolution  imaging  system,  the  problem  of  reducing  data 
handling  requirement  is  addressed  to  bit  compression  which  is  the 
minimization  of  the  number  of  bits  per  sample  used  to  represent 
the  image. 

Techniques  of  this  sort  are  usually  referred  to  a  more 
general  term  as  data  compression! 1 ] .  Unfortunately,  most  existing 
data  compression  techniques  are  found  not  to  be  useful  for  our 
application  in  microwave  imaging  because 

1)  In  applying  these  techniques,  both  the  input  and  the 
output  of  the  compression  processor  are  in  the  same  domain.  For 
microwave  imaging,  the  data  need  to  be  compressed  are  the  aperture 
data  which  are  the  complex  samples  of  the  electric  field  taken  by 
the  receiving  antenna,  which  is  usually  a  phased  array  (2).  The 
image  is  related  to  the  aperture  data  by  a  diffraction  integral  of 
physical  optics  13),  which  is  well  approximated  by  a  Fourier 
transform  when  the  distance  from  source  to  antenna  is  large  enough 
(4).  In  other  words,  the  expected  input  and  output  of  a 
compression  processor  are  in  different  domains  related  by  an 
integral  transformation.  It  is  the  nature  of  integral 
transformation  that  any  point  in  the  transform  domain  does  not 
correspond  to  a  point  in  the  source  domain,  but  to  the  entire 
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source;  As  a  result,  the  image  quality  criteria  which  govern  the 
data  compression  algorithms  are  no  longer  valid. 

2)  The  performance  of  data  compression  technique  is  measured 
by  the  capability  of  compression  (e.g.  compression  ratio  which  is 
the  ratio  of  the  number  of  bits  per  pixel  used  for  coding  the 
data-compressed  image  to  the  number  of  bits  per  pixel  used  for 
coding  the  original  image  ).  However,  from  image  fidelity  point 
of  view,  the  higher  the  compression  is,  the  larger  the  distortion 
is.  Therefore,  in  order  to  turn  these  techniques  into  a  useful 
application  in  image  processing,  one  requires  not  an  arbitrary 
high  compression  scheme  but  the  best  trade  off  between  the  image 
fidelity  and  the  degree  of  compression  for  given  system 
constraint.  In  other  words,  one  needs  to  know  how  the  image 
quality  is  related  to  the  controlled  parameters  of  the  compression 
processor,  so  that  for  given  image  quality  requirement,  one  can 
determine  what  minimum  number  of  bits  that  the  data  can  be 
compressed  to  or,  for  given  number  of  bits  limited  by  the  system, 
what  quality  the  constructed  image  can  achieve.  Unfortunately, 
such  theory  is  found  not  to  be  available  in  most  existing 
techniques . 

3)  Even  if  it  is  possible  to  develop  these  techniques  to  be 
adaptable  to  the  microwave  imaging  system  (meaning  that  the  first 
two  problems  are  resolved),  their  performances  are  still  not 
justified  in  considering  the  effect  of  low  correlated  source  of 
microwave  image  and  the  additional  cost  and  increasing  complexity 
of  hardware  and  software  in  applying  these  techniques. 

The  research  objective  was  to  develop  a  theory  of  Bit 
Compression  for  the  high  resolution  microwave  imaging  system  to 
determine  the  minimum  number  of  bits  required  of  the  quantized 
signal  for  the  image  to  achieve  any  given  quality  or,  conversely, 
what  quality  of  image  can  be  obtained  for  anv  given  number  of  bits 
that  the  quantized  signal  retains.  The  theory  will  be  developed 
by  1)  deriving  an  optimum  quantization  scheme  which  maximize  the 
image  quality  over,  the  entire  system  for  given  number  of  quantized 
level;  2)  deriving  a  theoretically  exact  and  physically  solvable 
analytic  relation  between  the  number  of  bits  per  sample  of  the 
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quantized  signal  and  the  quality  measure  of  image,  and 
analyzing  the  effects  due  to  data  compression. 
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2.  THE  THEORY 


2-1  The  System  Model 

A  typical  microwave  imaging  system  model  is  shown  in  Figure  1 

151  . 

Sir,1*1)  is  a  source  distribution  function  which  is  actually  a 
complex  reflectivity  function  due  to  illumination  from  a  distant 
microwave  transmitter.  The  electric  field  e(x)  induced  in  the 
aperture  is  related  to  S(r,$)  by  a  diffraction  integral.  An' angle 
$  is  measured  from  the  normal  to  the  aperture  of  the  linear  array. 
By  replacing  $  with  u  =  sin(i$i),  the  integral  transformation  ,  for  a 
far-field  source  or  a  near-field  source  with  a  well  focused 
antenna  array,  is  reduced  to  Fourier  transform  where  u  becomes  the 
canonically  conjugate  variable  to  x/x  in  the  Fourier  relation  [41. 

The  source  distribution  Str,#)  can  be  reduced  to  one 
dimensional  function  S^(u)  because  the  range  dimension  in  radar 
can  be  considered  as  a  discrete  set  of  range  bins  indexed  by  k 
1 5  ]  . 

w(x)  is  the  aperture  weighting  function.  The  current 
distribution  f(x)  in  the  aperture  is  the  product  of  e(x)  and  w(x). 
With  the  exception  of  modern  superresolution  techniques  [6],  the 
image  is  generally  the  magnitude  of  the  Fourier  transform  of  the 
current  distribution  151,  which  can  be  written  as  If  =  I  F(f)  I 
where  F(-)  denotes  Fourier  transform. 

When  the  data  compression  is  introduced,  the  compressed 
aperture  data  becomes  g  =  D  C  f  3  where  D  is  the  compression  operator. 
The  image  reconstructed  from  compressed  data  then  becomes  the 
magnitude  of  the  Fourier  transform  of  g  denoted  by  I  =  I F ( g ) I . 

The  selections  of  the  weighting  function  w(x)  will  not  change 
the  nature  of  image  distortion  due  to  data  compression  because  the 
reference  image  is  1^  in  which  the  effect  of  w(x)  has  been 
i ncluded . 

Finally,  the  problem  becomes  that  for  a  given  image  quality 
requirement  of  Ig,find  a  Bit  Compression  technique  to  compress  f, 
such  that  the  number  of  bits  of  g  can  be  minimized. 
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range 


Figure  1.  Microwave  Imaging  System  Model 


2.2  The  Signal  Model 

The  microwave  image  data  can  be  considered  as  a  random 
process.  By  definition,  a  random  process  is  an  ensemble  of  time 
functions  where  the  time  function  is  random  variable  for  a  given 
instant  of  time  173. 

The  random  process  model  for  microwave  imaging  is  illustrated 
in  Figure  2.  For  a  given  time  instant  which  corresponds  to  one 
snap  shot,  the  image  data  is  the  spatial  sample  of  the  radiation 
source  for  one  range  bin,  therefore  the  time  functions  are  the 
receiving  signals  in  the  aperture  and  the  ensemble  is  defined  to 
include  all  array  elements. 

The  process  is  assumed  to  be  ergodlc  because  images 
corresponding  to  different  snap  shot  are  assumed  to  be 
uncorrelated . 

The  process  is  not  stationary  in  general.  However,  the 
statistical  study  of  data  distributions  shows  that  if  scenes  to  be 
imaged  are  divided  into  different  groups  according  to  different 
categories  of  objects  (e.g.  residence,  industrial, 
airplane, ... ,etc)  then  the  process  can  be  considered  as  stationary 
within  a  group  provided  that  their  data  distributions  are  close 
enough  and  the  parameters  of  the  processor  are  not  sensitive  to 
variations  of  the  data  distributions  in  the  group. 

Therefore,  it  will  be  assumed  that  both  the  stationarlty  and 
the  ergodicity  are  held  for  the  model. 

2.3  Basic  Definitions 

Image  Fidelity  Measure:  Define  the  normalized  zero-lag 
correlation  coefficient  p  between  the  original  and  the  distorted 
image  as  a  measure  of  image  quality: 


II  - 

frms  grms 

where  I^denotes  the  original  image  and  I^denotes  the  image 
obtained  from  the  distorted  ( i .e . , quantized )  input  data.  Both 
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I,and  Iare  functions  of  the  image  coordinates.  The  overbar  means 

f  g 

average  over  all  picture  elements  in  the  images  and  rms  denotes 
root  mean  square. 

In  most  image  processing  applications,  the  mean  square  error 
e  is  used  to  measure  the  image  fidelity  [1].  It  has  been  shown 
(51  that  s  and  p  are  equivalent  if  s  is  normalized  and  if  the 
distorted  signal  g  is  scaled  such  that  its  energy  and  the  energy 
of  the  original  signal  f  are  equal.  However,  these  two  major 
differences  do  make  the  correlation  coefficient  metric  more  useful 
in  our  application  than  the  conventional  mean  square  error  metric. 
Because  of  the  scaling  property,  the  comparison  between  the  proper 
image  and  the  distorted  image  is  made  under  the  same  energy  base, 
so  that  only  the  unsimilarity  will  be  detected  from  the 
measurement.  Because  of  the  normalization  property,  the 
sensitivity  of  the  image  quality  measure  to  the  variation  of  data 
distribution  for  different  scene  will  be  significantly  reduced, 
and  a  bounded  relative  error  scale  will  usually  make  more  sense 
than  an  unbounded  absolute  error  scale  e  . 

The  selection  of  correlation  coefficient  metric  become 
significant  in  the  optimum  quantization  procedure  where  it  leads 
to  Independent  quantizer  design  while  the  mean  square  error  metric 
lead  to  dependence  between  amplitude  quantizer  and  phase 
quantizer . 

Information  Measure:  Define  entropy  H  to  measure  the 

information  content  of  the  quantized  signals; 

HZ  ’  PZl  (2'2) 

where  Z=A  representing  amplitude,  Z=£  representing  phase  and  pzils 
the  probability  that  the  signal  z  takes  a  value  in  the  ith 

quantization  region  . 

2.4  Aperture  Domain  Formulation 

We  assume  that  the  source  or  scene  being  imaged  is  in  the 

far  field.  p  can  then  be  expressed  in  the  aperture  domain  by 
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using  the  Fourier  transform  relationship  between  the  source 

distribution  and  the  electric  field  induced  in  the  aperture,  and 

Patseval’s  Theorem  equating  energies  in  the  two  Fourier  domains 

[5).  In  the  following  derivation,  f  consists  of  the  measured 

samples  of  the  electric  field  in  the  aperture  and  g  are  its  values 

after  quantization.  f  and  g  are  functions  of  the  aperture 

coordinates.  I*«|F(f)|  and  I  =|F{g)|.  F(-)  denotes  Fourier 

f  g 

Transform  and  *  denotes  the  complex  conjugate. 

From  (2.1), 


P  s 


* frms1 grms 


Because 


|F( f ) |  |F(g) 1 


[  F(f)  1rBJ  F(g>  1 


rms 


we  write 


I  Ft  f ) I  I F ( g ) |  *  | Ft  f )  F(g) 


I  Ft  f ) I  | F ( g ) |  =  K 


F(f)  F ( g ) 


where  K  greater  than  or  equal  to  1.  Then 


p  =  K 


F  (  f )  F  ( g ) 


(F(f) J  (F(g)  1 

rms  rms 


=  K 


F ( f )  F(g) 

C  I F ( f  )  ?  I F(g)  f  ]l/S 


Parseval's  theorem  implies  that 
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F(f)  F(g) 


C  I  F(  f )  f  IF(g)  f  ]1/S  Cl  f  f  I  g  f]1/2 


Therefore 


p  =  K 


& 

=  K  P, 


(2.3) 


£  g 
rms’rnis 


The  analysis  is  performed  in  polar  coordinates.  Let  f=Ae 

and  g=A'e  .  Reference  [5]  shows  phase  and  amplitude  to  be 
independent  when  the  scene  being  imaged  consists  of  many  sources 
and  none  of  the  sources  dominates  all  the  rest.  Assuming  this 
condition,  PQ  can  be  expressed  as 


po  = 


A  A 1 
rms  rms 


where 


s  p  P 

A 


PA  =  A  A'/  A  A' 

A  rms  rms 


exp  <j(£-S')) 


(2.4) 


(2.4a) 


(2.4b) 


where  the  overbar  means  ensemble  average  over  all  array  element 
which  could  be  replaced  by  the  expectation  value  because  the 
number  of  array  elements  is  very  large. 

2.5  Optimum  Quantizer 

Techniques  of  optimum  quantization  usually  deal  with  real 
input  sequences.  Examples  are  the  early  papers  by  Max  18]  and 
Paez  and  Glisson(9].  In  high  resolution  microwave  imaging,  the 
aperture  data  to  be  quantized  are  complex  random  variables.  In 
this  thesis,  an  optimum  quantization  scheme  for  complex  data  in 
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polar  coordinates  Is  introduced  .  Gallagher  (10]  also  dealt  with 
optimum  quantization  of  complex  data  in  polar  coordinates.  There 
are  two  main  differences  between  these  two  papers  and  this  work. 

First,  in  (81,(9]  and  (10],  only  the  quantizer  was  optimized, 
whereas  in  our  procedure,  the  optimization  is  carried  out  over  the 
entire  imagine  system  (which  includes  the  transform  from  source 
distribution  to  electric  field  in  the  aperture  and  the  inverse 
transform  from  quantizer  output  to  image). 

Second,  the  optimization  procedure  in  [10]  leads  'to  an 
amplitude  quantization  scheme  that  is  dependent  upon  the  phase 
quantization.  In  our  procedure,  which  maximizes  the  image 
correlation  coefficient  P  (2.1),  the  amplitude  and  phase  quantizer 
are  independent  This  is  significant  because  1)  it  gives  the 

complete  freedom  to  achieve  the  most  efficent  bit  allocation 
between  amplitude  quantizer  and  phase  quantizer;  2)  it  makes  it 
possible  to  express  the  total  entropy  of  the  quantized  signal  with 
the  individual  entropy  measurement  from  the  two  quantizers'  output 
without  involving  conditional  entropy  [7]  which  is  usually 
difficult  to  evaluate;  and  3)  it  leads  to  a  similar  optimum 
conditions  both  for  amplitude  quantizer  and  phase  quantizer  as 
those  developed  for  real  signal  by  Max,  so  that  one  is  able  to  use 
all  resources  and  facilities  developed  for  Max  quantizer,  such  as 
Bruce's  algorithm  111]  and  Sharma's  algorithm  (12],  to  calculate 
the  optimum  settings. 

The  conditions  of  optimum  quantization  will  be  obtained  by 
maximizing  the  image  correlation  coefficient  P. 

It  can  be  found  later  that  the  conditions  which  maximize  Pn 
also  maximize  K.  Therefore,  the  image  correlation  coefficient  P 
will  be  maximized  as  the  lower  bound  correlation  coefficient  PQ  is 
maximized . 

Because  p,-,  =  p.  p,  ,  and  because  p.  and  pt  are  independent, 
U  A  •£  A  i 

pQ  is  maximized  as  each  of  pa  and  P^  is  maximized  individually. 

2.5.1  Optimum  Amplitude  Quantizer 

Let  be  the  number  of  quantized  levels,  aj  be  the  ith 

output  level  of  the  quantizer,  the  value  assigned  to  a  sample 


11 


falling  in  the  interval  [  ajL_1/al  ),  and  be  the  endpoint  of 
input  range  of  the  quantizer  with  i=l, 2, . . . ,NA(z  replaced  by  a  in 
Figure  3). 

a 

Define  a^=  E(  A  I  a^_^  <  A  <  a^)  with  E(-)  representing 
expectation  value,  pM*  Pr(  a^  <  A  <  ai )  with  Pr(-)  standing 
for  probability. 

The  terms  in  (2.4a)  may  be  calculated: 


N. 


A  A1 


*  2  Prta^  <  A  <3i)  E  (  AA'Ia^  <  A  <  a£ ) 


and 


NA 

=  S  PMa;  E(  *1  aj.j  <  »  <  aj) 


^  A 

*  p*i»i  n 


Arms=  C  ^  pA1E(A'  1  al-l  <  A  <  ai)V 


\  2l/s 

’  C  fpMai  5 


The  amplitude  correlation  coefficient  is 


NA 

\  PAlalal 


NA 

C  pAiai 


,  V2 

“  3  A 


rms 


(2.5) 


It  can  be  shown  that  is  maximized  if  aj  and  a^  satisfy 


a-=  at 


i=l,2, . . .,N, 


(2.6a) 


ai  = 


-V  ai+i 


i=l, 2, . . . N. -1 


(2.6b) 
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'igure  3.  Description  of  Definitions  of 
Quantizer  Design 


2.5.2  Optimum  Phase  Quantizer 

Let  be  the  number  of  quantized  levels,  4J  be  the  output 

level  of  the  quantizer,  the  value  assigned  to  a  sample  falling 
within  the  interval  **i-i'*i^#  and  *i  be  the  en<3Point  of  the 
input  range  of  the  quantizer  with  i=l,2, . . .  ,N^-l(z  replaced  by  4 
in  Figure  3) . 

A 

Define  4^  =  E(  5  I  4^_^  <  2  <  4^  )  E(-)  representing 
expectation  value,  and  p^j  Pr(  4^_^  i  2  <  4^  )  with  Pr(-) 
standing  for  probability. 

The  phase  correlation  coefficient  (2.4b)  becomes 

p£=  |  expl j (2-2* ) 1  I 

=  Is5  p51E{  exp( j(2-2' ) ] I  4i_1<  2  <  4 ± ) I 

It  can  be  shown,  for  uniformly  distributed  phase,  that  is 
maximized  if  4^  and  4J  satisfy: 

4 |  1=1,2, . .  .,Ng  (2.6c) 

A 

^i+  ^i+1 

4,  =  — 1 - i-1,2,  .  .  .,N,-1  (  2 . 6d  ) 

1  2  * 

Equations  ( 2 . 6a ) - ( 2 . 6d )  are  the  design  equations  of  the 

optimum  quantizer.  From  (2.6),  it  is  clear  that  the  dependence 
between  the  amplitude  quantizer  and  the  phase  quantizer  in 
Gallagher's  design  vanishes.  The  equations  ( 2 . 6a) -(2 .6d)  say  that 
in  the  optimum  design  of  the  two  quantizers,  the  ith  output  level 
equals  the  conditional  mean  of  the  input  data  in  the  ith  range  and 
the  endpoint  of  the  ith  range  equals  the  arithmetic  mean  of  two 
adjacent  output  levels.  These  statements  are  true  for  arbitrary 
amplitude  distribution,  for  uniform  phase  distribution  and  for 
amplitude  and  phase  being  independent. 
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2.6  The  Transfer  Index 

From  (2.1)/  one  can  see  that  PQ  has  the  same  form  as  P  except 
that  It  Is  defined  as  the  magnitude  of  complex  quantity. 
Therefore  we  can  consider  K  as  a  transfer  index  which  transfers 
the  image  correlation  coefficient  from  the  aperture  domain  into 
the  image  domain.  The  expression  for  K  can  be  derived  in  the 
following  procedure. 

From  (2.1)/ 


p  = 


I-  I 
f  9 

I.  I 
frms  grms 


I F ( f ) |  I F(g )  I 

(  F(f)  )  I  F(g)  ] 

rms  rms 


But  from  (2.3) 


pO= 


F (  f )  F(g) 


f  g 
rms^rms 


(F(f) ItF(g) 1 

rms  rms 


where 


=C  I  |F(f ) |exp( j*f ) |F(g) |exp(-j*g) 

C  =  {  (  F(f)  1  (  F(g )  ll"1 

rms  rms 


Independence  between  phase  and  amplitude  In  the  image  domain  (5) 
implies  that 


P0=  C  I F ( f ) |  ! F ( g ) |  lexp( j(*f-*g) ) 


=  lexp( 3 (+f-*) ) I  P 
which/by  definition,  gives 


K  =  I  exp(  j  (  Pf-Py  )  ) 
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Let  A^  denotes  the  measured  amplitude  in  the  ith  element,  $ ^ 
denotes  the  measured  phase  in  the  ith  element,  and  n  '  M 
represents  values  after  distortion.  After  much  mathematics,  it  is 
found  that 


-l 


K  =  exp(jcc) 


where 


a  =  tan  *6 


and 


(Z  A.sin4\)(Z  A!cosT')-(Z  A,cost.)(Z  A'sint'  ) 


|3= 


A 


A 


A 


A 


L 


_L 


X 


_L 


(Z  A1cos,f1)(Z  A'cost|)  +  (Z  A1slnfi)(Z  A'sintp 


(2.7) 


where  with  v^being  the  Fourier  kernel,  and  the 

summation  is  taken  over  all  elements  (i  from  1  to  N). 


*  .-1? 


By  defining  V  =  N  Z  V^,  e  can  be  expressed  as 

i 


(  A  sinf  )(  A'costM-t  A  cos')'  )(  A'sint') 


(2.8) 


(  A  cosT  )(  A'cos'V,)  +  (  A  sint  )(  A' si  n't”) 


Because  the  number  of  elements  is  very  large,  V  can  be 
approximated  by  the  expectation  value  of  X  .  This  Implies  that 
the  transfer  index  K  is  uniquely  specified  for  given  distribution 
of  aperture  data  and  the  quantizer  parameters  as  the  number  of 
elements  goes  to  very  large. 

This  property  shows  that  the  exact  solution  of  image 
correlation  coefficient  exist  if  the  distribution  of  aperture  data 
are  given.  It  also  allow  us  to  obtain  K  by  using  computer- 
simulated  Images  through  the  relation  K=p/pq  with  p  obtained  from 
(2.1)  and  PQ  obtained  from  (2.4).  This  method  is  more  accurate 
than  calculating  (2.7)  because  of  the  tremendously  reducing  number 
of  mathematical  procedures  that  are  involved;  the  result  is 
smaller  accumulation  of  errors  in  the  calculation. 


In  section  2.5,  it  was  shown  that  PQis  maximized  when  the 
optimum  quantization  conditions  are  applied.  A  numerical  study  of 
(2.7)  shows  the  same  result  for  K,  namely  that  K  is  maximized  when 
the  optimum  quantization  conditions  are  applied.  Furthermore  it 
was  found  that  K  is  insensitive  to  changes  in  the  output  levels 
provided  that  the  maximum  deviation  from  the  optimum  set  is  less 
than  10  per  cent. 

Some  values  of  the  estimate  of  K  are  listed  in  Table  1. 

A  numerical  study  of  K  based  on  (2.7),  and  the  results  of 
experimental  testing  discloses  the  following  additional  important 
properties  of  K: 

1)  Using  the  estimates  of  K  into  the  evaluation  of  image 
correlation  coefficient  with  experimental  data  (Tables  5  and  6) 
infers  that  K  is  insensitive  to  variations  in  the  distribution  of 
the  complex  aperture  data.  This  property  is  very  important  in  a 
practical  sense  because  it  indicates  that  K  is  relatively  scene 
free. 

2)  The  transfer  index  K  is  more  sensitive  to  phase  distortion 
than  to  amplitude  distortion.  As  Table  1  shows,  K  is  closer  to  1 
for  amplitude-distortion-only  than  for  phase  distortion. 

3)  The  same  table  shows  that  when  the  phase  distortion  is 
significant  (e.g.,  1-2  bits  of  phase),  K  is  determined  almost 
entirely  by  the  phase  quantizer. 

4)  When  the  phase  distortion  is  insignificant  (e.g.,  3  or 
more  bits  of  phase),  K  is  determined  almost  entirely  by  the 
amplitude  quantizer. 

2.7  The  Analytic  Equations 

By  applying  the  optimum  condition  (2.6)  to  (2.4)  and 
summarizing  the  results  from  previous  sections  we  obtain  the 
analytic  equations  of  image  correlation: 


H  (bits) 


Table  1 .  Estimate  of  K  .  The  ranges  R1  and  R 2  defined  by 

<  .  <  <  .  < 

R1 : (0. 6  •  pa  •  0.9)  and  R2:(0.9  •  pA  •  1.0)  where  pa-  A/Arms 


(2.9a) 


P.  =  a_m_ (a. ) /A 
A  r ms  i  rms 


N, 


1/2 


Vhere  aris*  C  ^  PAlai  ] 


'H$ 

p,=  sinc(n  2  *) 

a 


(2.9b) 


NZ  '  ‘ 

HZ(zi)  =  ‘  l  PZi(2i)  log  Pz^V 

vhere  Z=A  representing  amplitude  and  Z =$  representing  phase. 

The  independence  between  the  amplitude  quantizer  and  phase 
quantizer  implies 


H  =  H.  +  Ht 
a  * 


K  is  determined  by 


K  =  exp(jco) 


-1 


vhere 

and 


cc  *  tan  13 


(  A  sint  )(  A'cos't")-(  A  cos'f  )(  A'sint') 


tS  = 


(  A  cos'V  )(  A'cosT')+(  A  sin1)1  )(  A'sinT ') 


For  given  system  constraints,  this  set  of  equations  relates 
the  optimum  image  quality  to  the  minimum  number  of  bits  per  sample 
of  the  quantized  signal  through  the  optimum  endpoints  setting  (z.) 

a  A 

and  optimum  output  levels  {z^} 

Table  2  lists  some  results  of  exact  image  correlation 
coefficient  P  versus  the  number  of  bits  per  sample  in  amplitude 
and  phase  for  Rayleigh  distributed  amplitude  and  uniformly 
distributed  phase  which  shows  that  the  image  quality  is  well 
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r°" 

] 

H^(bita) 

2 

~i 

3  1 

i 

ha 

0.806 

0.886 

1 

0.9*19  | 

(bits) 

i  i 

i 

0.873 

0.961 

0.985  j  p-Kp. 

i 

i  2 

i 

_i _ 

0.896 

0.986 

0.990  | 

1 

_ i 

Table  2.  Image  Correlation  Coefficient  vs.  the  number  of 
Bits  Per  Sample:  Amplitude  and  Phase 


preserved  for  high  compressed  aperture  data  (for  example,  from  the 
first  column  first  row  of  Table  2,  one  can  see  that  more  than  80% 
of  the  image  fidelity  is  achieved  for  only  1  bit  of  aperture 
data ) . 

As  one  can  observe  that  both  the  analytic  equations  of  image 
quality  and  the  optimum  conditions  have  no  closed-form  solution. 
They  can  only  be  solved  numerically  by  using  large  computer 
program.  However,  in  order  to  study  the  properties  of  the  theory, 
it  is  not  always  necessary  to  carry  out  the  general  and' exact 
solutions.  Good  approximation  and  simulation  using  typical 
distribution  are  usually  more  desirable. 

2.8  Suboptimum  Conditions 

It  is  well  known  that  phase  carries  more  information  useful 
for  image  construction  than  amplitude.  Therefore,  a  suboptimum 
amplitude  quantization  schemes  could  yield  near-optimum  results 
provided  the  phase  is  optimally  quantized. 

If  one  replaces  the  condition  of  (2.$b)  with  the  condition  of 
setting  all  pAi  to  be  equal: 

'ha 

PAi=  p  =  2  *  for  all  i  (2.10) 

then  the  quantization  procedure  can  be  reduced  to  one  dimensional 
and  P'A  can  then  be  expressed  explicitly  as  the  function  of  HA 

HA 

-H  2  1/2 

P.  »  C  2  2  CT ( H _ )  ]  (2.11) 

A  i=l  1  A 


where 


a \/  Arms 


Because  C^is  only  a  function  of  HA  ,  pa  is  uniquely  specified 
by  Hk. 
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It  is  found  by  computer  calculation  of  (2.9a)  and  (2.10)  that 
P^(2.9a)  is  almost  identical  to  p^(2.10)  in  the  sense  that  the 
largest  error  between  them,  which  appears  when  NA=2,  is  less  than 
0.2*.  However,  equation  (2.10)  greatly  simplifies  the  analysis 
procedure . 


2.9  Simulation 

Assuming  that  the  scene  to  be  imaged  is  in  the  far  field  and 
consist  of  many  sources  (targets)  characterized  by  a'  real 
amplitude  A^  and  phase  ,  the  electric  field  Induced  in  the 

antenna  array  is  e(x)  =  2  A^  exp(  jkxu^  +  6^  ).  Further  assuming 
that  none  of  the  sources  dominates  the  sum  of  all  the  rest  which 
suggests  that  all  these  sources  are  Independently  identically 
distributed.  By  the  Central  Limit  Theorem,  the  quadrature 
components  of  the  electric  field  are  independent  Gaussian 
processes  which  implies  that  the  amplitude  and  the  phase  of  the 
current  distribution  f  are  Independent  with  Rayleigh  distributed 
amplitude  and  uniform  distributed  phase.  By  using  this  multiple 
targets  model  (51  and  applying  the  suboptimum  condition  (2.10),  a 
closed  form  solution  of  lower  bound  image  correlation  coefficient 
PQln  terms  of  the  entropy  of  quantized  signal  can  be  obtained. 

To  see  this,  recall  equation  (2.11) 


-H 


"V  <  2 


A  1/2 

2  Cj  )  ' 


where 


with 


C,  =  a. /A 
i  i  rms 


-H, 


e*r  t  "2 


for  all  1. 


By  Noting  that  the  pdf  of  amplitude  is  Rayleigh,  that  is 


PA(a)  =(a/a“)exp(-(a"/2a“) ) 

Evaluating  the  terms  in  (2.11)  gives 

at 

A„  -  (  X  a2P.(a)  da  )/Z  =  A  a  (2.12) 

rms  q  A 

where  o  is  the  standard  deviation  of  amplitude. 
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E(A|a^_^<  A  <  a^  )  =  p  l<f  a  P^a)  da 

ai-l 

*  p  1C-a1exp(-a^/2a2)+a1_1exp(-a^_1/2CT2) 

+  V-2TT  a  [Erf  (at/CT)-Erf  (a^j^/cr) }} 

where  Erf(x)  is  the  error  function  of  x  defined  by 

x 

Erf (x)  =  {v'2tt  )  1  J  exp(-y“/2)  dy 

—to 

After  much  mathematics,  a^s  found  to  be: 

al*{2  ln[  (l-ip)"1]} l/a  1*0,1, ... ,NX-1 

and 

aM*  >» 

NA 

Now  let 

CjMj/Kn/2)  l/a] 
and  note  that 


then 

Cj*  2  {  InE  ( 1-12  A)_1  }/tt}  1/Si«0,l,  .  . -,Na-1 

Substitute  (2.15)  into  (2.13)  and  then  substitute  the 
(2.13)  and  (2.12)  into  Ct  in  (2.11).  The  result  is 

HA 

Ct=  2  ,tt  {  0.5  (  C1_1exp(-nC^_1/4)-Ciexp(-TTC^ 
+  C  Erf  (C^n/2)  -  Erf  (Ci_1(H/2)  1/2)  ]} 


(2.13) 


(2.14) 

(2.15) 

(2.16) 

mod i f  ied 

/A)  ) 
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(2.17) 


where  1»1# 2, . . . ,NA~1  ,  and  C^is  determined  by  (2.16). 
A  summary  of  results  follows: 


where 


H  =  H,+  H, 

A  S 

HA 

-H  2  1/2 

P.  -  C  2  A  S  C“(H.)  ] 
a  i=l  1  A 

HA  -  2 

Ci«  2  vTT  {  0.5  (  C1_1exp(-nCi_1/4)-Ciexp(-nCi  /A)) 

+  [  Erf  (Ci(n/2)  Vf  -  Erf(C1_1(n/2)  1/2)  ]} 

•  || 

C1=  2  {  In [( 1-12  A)“1]/TT)1/2i*0/l#.../NA-l 


pa  can  be  approximately  given  by 


Lastly 


p  =  l-exp[-(1.056H?,95+2.173) 1 

A  A 

-H* 

P,=  slnc(n  2  *) 


Results  from  carrying  out  these  equations  are  shown  In  Table 
3,  Table  A,  Figure  A  and  Figure  5.  Table  3  lists  pq  for  various 
combinations  of  H.  and  HT.  Table  4  relates  the  correlation 

A  '5 

coefficient  PQ  to  a  subjective  assessment  of  image  quality  [5]. 
Figure  4  plots  PAand  against  the  number  of  amplitude  and  phase 
bits  per  sample  and  Figure  5  plots  pq  where  the  dashed  lines  are 
for  constant  sums  of  H,  and  H_. 

A  'jL 

An  Important  observation  drawn  from  the  analytic  equations 
( ( 2 . 11 ) - ( 2 . 17 ) )  of  lower  bound  correlation  coefficient  Pn  is  that 
pQ  is  independent  with  the  statistical  properties  (mean  or 
variance)  of  the  image  data  and  uniquely  related  to  the  entropies 
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p. 

Image 

0.6 

often  recognizable 

0.8 

often  acceptable 

C  .9 

good  quality 

0.99 

high  fidelity 

Table  Experimental  relation  of  correlation  to  linage  quality 


1 


2 


4 


5 


re  4* 


3 

H  AND  H,  (3ITS) 
A  9 


vs.  Number  of  Amplitude  Bits  Per  Sample  and 
vs.  Number  of  Phase  Bits  Per  Sample 


2.7 


Figure  5.  Lower  Bound  Correlation  Coefficient  As  a 
Function  of  The  Number  of  Amplitude  Bit  HA 

and  Phase  Bits  H.  Per  Sample 


of  the  quantized  signals.  This  conclusion  can  also  be  considered 
as  correct  for  the  exact  inage  correlation  coefficient  p  because 
the  transfer  index  K,  when  the  number  of  phase  bits  is  small,  is 
asymptotically  independent  with  amplitude  which  is  the  one  being 
the  function  of  mean  and  variance. 

It  is  shown  that  the  phase  carrying  more  information  useful 
for  image  construction  than  amplitude.  This  can  be  seen  in  Figure 
4  which  shows  that  dominates  the  value  of  pq. 

It  is  also  shown  that  hard-limiting  the  input  signals  to  a 
phased  array  while  leaving  the  phase  unaltered  preserves  much  of 
the  image  integrity.  This  agrees  with  work  by  Kermish  [15),  Van 
Hove  [16]  [171,  Oppenheim  [18),  and  Steinberg  [5).  The  result  is 
seen  in  Figure  5  on  the  upper  portion  of  the  H^axis  (HA=0):  the 
correlation  coefficient  PQ  in  this  region  is  nearly  0.9  indicating 
(Table  4)  good  image  quality.  Amplitude-only  inf ormation, on  the 
other  hand,  destroys  the  image.  Hard  limiting  plus  one  bit  of 
phase  information  sometimes  suffices  (PQ*  0.6).  Hard  limiting 
plus  two  bits  of  signal  phase  per  sample  results  in  a  correlation 
coefficient  of  0.8,  vMch  often  results  in  acceptable  imagery. 

Figure  5  shows  that  a  minimum  of  four  bits  of  information  per 
complex  signal  sample  guarantees  a  relatively  good  image  (  SO. 9) 
provided  that  three  bits  are  used  for  the  phase  and  one  bit  for 
the  amplitude.  Eight  bits  are  required  to  ensure  a  high  fidelity 
image  (  SO . 99 ) . 

Table  3  and  Figure  5  provide  a  quantitative  rule  of  thumb  for 
assigning  data  bits  to  phase  and  amplitude:  When  the  number  of 
bits  per  data  sample,  H.+  Hi;  is  no  more  than  8,  the  number  of 

A  be! 

bits  assigned  to  phase  should  be  two  larger  than  the  number  of 
bits  assigned  to  amplitude.  When  H.+  H_>.  10,  H  =  H^is  the  optimum 

A  £  A  !*! 

assignment.  These  rules  also  show  that  the  phase  information 
thoroughly  dominates  the  process  which  is  again  agree  with  the 
well  known  experimental  results  in  coherent  imaging  (e.g. 
holography  or  microwave  imaging). 


3.  EXPERIMENTS 


We  have  conducted  6  experiments  to  test  the  theory.  The 
experimental  data  are  obtained  from  near-field  tvo-dimensional 
high-resolution  radio  camera  images  having  3m  of  range  resolution 
and  2-6m  cross-range  resolution  C2J.  The  wavelength  is  3cm.  The 
experiment  numbers  and  the  target  descriptions  are  as  follows: 

301  Limerick  Nuclear  Power  Plant 

303  Cromby  Power  Plant 

308  Marshall  Street,  Residential  Houses 

309  High  Street,  Residential  Houses 

310  Potato-Chip  Factory 

312  Suburban  Housing  Development 

The  targets  are  all  in  Phoenlxville,PA,  at  distances  of  4.5 
to  17  km  from  the  radio  camera,  which  is  located  at  the  Valley 
Forge  Research  Center  of  the  University  of  Pennsylvania.  The 
image  correlation  coefficient  obtained  from  the  experimental  data 
is  denoted  by  px  which  is  obtained  by  evaluation  of  (2.1)  using 
images  produced  by  experimental  data  with  optimum  quantization 
applied. 

The  value  predicted  by  the  theory  is  p=KpCj  where  pq  obtained 
by  (2.9a)  and  (2.9b)  using  measured  data  in  the  aperture  and  K  is 
obtained  from  Table  1. 

3.1  Phase  Distortion  Only  Experiments 

In  this  case  p=Kpt  because  P  =1.  pt  was  directly  calculated 

a  A  j t 

from  (2.9b).  The  comparison  between  Px  and  P  for  the  six 
experiments  are  shown  in  Table  5.  The  seventh  row  shows  the 
average  of  the  measured  correlation  coefficients  of  the 
experimental  data.  In  the  eighth  row  are  the  predicted  values. 

From  Table  5  one  can  observe  that  the  theoretical  value  is 
very  close  to  the  average  experimental  value.  This  is  a  highly 
practical  result  because  the  theoretical  correlation  (2.9b)  was 
derived  under  the  assumption  of  uniformly  distributed  phase.  In 


30 


run 

r 

i 

I 

1  bit  of  phase  * 

2  bits  of  phase 

Remarks 

301 

0.89  j 

0.98 

303 

0.92  1 
| 

0.93 

308 

0.91  j 

0.98 

Experimental 

Px  309 

0.90  ? 

0.98 

correlations 

310 

0.92  j 

0.93 

312 

0.93  1 

0.98 

1 

0.91  1 

0.98 

Ave.  exp. cor re 1. 

P-Kp,-Kp^ 

0.91  | 

1 

1 

_ L 

0.99 

_ 1 

Theoretical  value 

Table  5  p^  Coopered  with  p.  Phase  Distortion  Only 


Axiauth  (HHlir«41«ni) 


actuality,  the  real  phase  distributions,  obtained  by  study  of  the 
phase  histograms,  are  not  quite  uniform  and  are  d 1st inguishably 
different  from  scene  to  scene.  This  means  that  P  Is  Insensitive 
to  variations  in  the  distribution  of  phase. 

The  value  of  P  for  the  *no  phase*  case  turns  out  to  be 
Indeterminate  because  (2.8),  from  which  K  Is  determined,  become 
Indeterminate.  It  was  shown  experimentally  In  (5)  that  no-phase 
data  produce  a  meaningless  Image  (Figure  6).  It  Is  for  this 
reason  that  this  case  is  ignored  in  Table  5. 


3.2  Amplitude  Distortion  Only  Experiments 

In  this  case  P=KP.  because  P_=l.  The  theoretical  values  are 

A  2 

given  in  the  left  column  of  each  group  of  data  in  Table  6.  K  was 
obtained  from  Table  1.  P.was  calculated  from  (3.3)  in  which 

A  A 

A  and  the  a,  were  estimated  from  the  individual  data  sets, 
rms  i 

The  second  column  is  P^ . 

The  third  column  is  P^  ,  which  is  obtained  from  images 
produced  with  the  optimum  quantizer  design  based  on  the  assumption 
that  each  distribution  was  Rayleigh. 

Four  observations  can  be  drawn  form  Table  S: 

1)  The  values  of  P  predicted  by  the  theory  are  very  close  to 
Px  obtained  from  the  images. 

2)  Unlike  what  was  observed  in  the  phase-distortion-only 
experiment,  the  image  correlations  are  remarkable  different  from 
scene  to  scene  specially  when  the  number  of  amplitude  bits  is 
small  (  e.g.,  for  hardlimiting,  P..,ni  =  0.69  while  P.7OS*0 . 89,  and  all 
of  them  are  significantly  smaller  than  the  image  correlation 
obtained  from  Rayleigh  distributed  amplitude  data,  which  is 
p=0.94).  One  can  also  see  that  image  correlation  coefficients  are 
similar  when  their  amplitude  histograms  are  similar,  (e.g., data 
sets  of  runs  309  and  310  and  of  308  and  312.)  If  we  define  the 
image  point  occupancy  to  indicate  the  "mass"  of  the  image: 
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o 


/  I 


f  rms 


(3.1) 


33 


1 

o  ! 

HA(bits) 

r 

i 

i  i 

_ i_ 

2 

run 

mam 

Emm 

maam 

301 

0.70  0.70  0.69  1 

0.80  0.81  0.84  i 

0.88  0.90  0.89 

303 

| 

0.79  0.77  0.77  | 

0.92  0.91  0.92  | 

0.97  0.97  0.96 

308 

0.88  0.89  0.89  1 

j 

0.95  0.95  0.95  * 

| 

0.96  0.97  0.97 

309 

0.78  0.79  0.79  | 

0.89  0.89  0.91  | 

0.93  0.95  0.95 

310 

0.30  0.81  0.81  | 

0.92  0.90  0.91  1 
| 

0.95  0.95  0.94 

312 

0.85  0.86  0.86  | 

0.93  0.94  0.94  | 

0.96  0.97  0.96 

Rayleigh 

0.95  0.94  0.94  ' 

1 

1 

0.99  0.98  0.98  | 

1 

0.99  0.99  0.99 

Table  6.  Px  Compared  with  p.  Amplitude  Distortion  Only,  p  is  the 
image  correlation  predicted  by  the  theory.  Px&nd  p£  are  the  measured 
correlations  when  the  quantizer  design  was  optimized  for  each  measured 
aperture  distribution,  and  for  the  Rayleigh  distribution,  respectively. 
The  values  in  the  bottom  row  are  based  on  simulation  data  having  a 
Rayleigh  distribution. 
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which  is  the  normalized  mean  value  of  image,  or  more  imaginally 

the  image  correlation  between  the  image  and  a  homogeneous  black 

image  (hardlimiting  image).  The  results  of  P  versus  I  are 

tabulated  in  Table  7.  From  Table  7,  one  can  see  that  P  increases 

toward  the  theoretical  value  as  I  increases. 

o 

3)  It  is  evident  that  the  estimates  of  K  from  Table  1  which 
are  used  to  evaluate  p  in  Table  6  work  quite  well  for  all  scenes 
although  all  of  them  have  different  amplitude  distributions.  This 
means  that  K  is  nearly  a  scene-free  parameter. 

4)  A  very  important  observation  from  the  comparison  of  Px  and 
p^  is  that  the  quantizer  design  based  on  the  Rayleigh  distribution 
gives  nearly  identical  performance  ,  in  almost  all  cases,  to  the 
designs  using  the  six  actual  distributions.  This  implies  that  one 
quantizer  design  is  appropriate  for  all  scenes. 


run  # 


^(Hardlimiting) 


301 

303 

309 

310 
312 
308 


Table  7. 
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I 


| 

i 

t 

! 

1 

i 


0.70 

0.37 

0.77 

0.M2 

0.79 

O.H2 

0.81 

0.50 

0.86 

0.55 

0.89 

0.58 

P^ts.  Inage  Points  Occupancy 


4.  DISCUSSION  AND  CONCLUSION 


4.1  Practical  Considerations 

Observations  drawn  from  the  experimental  testing  o£  the 
theory  can  be  summarized  into  the  following  categories  with 
practical  considerations,  and  two  major  aspects  of  the  theory, 
which  are  the  optimum  quantizer  design  and  the  estimate  of  image 
correlation  coefficient,  will  be  examined. 

1) Sensitivity:  The  qualitative  analysis  and  numerical  study 
shows  that  the  transfer  index  K  is  insensitive  to  the  amplitude 
distribution  and  to  the  variation  in  phase  distribution,  while 
Table  6  shows  that  PQ  is  relatively  sensitive  to  the  changes  of 
amplitude  distribution.  The  optimum  quantizer  design  is  found  not 
sensitive  to  the  data  distributions. 

2) Generality:  The  theory  was  developed  under  the  assumption 
of  a  uniform  distribution  of  phase  in  the  aperture  data;  no 
restriction  was  Imposed  upon  the  amplitude  distribution.  However, 
since  the  variation  from  uniform  distribution  that  the  phase 
distributions  of  actual  data  have  is  considered  as  small  and  both 
optimum  quantization  conditions  and  estimate  of  P  are  insensitive 
to  the  phase  variations.  The  theory  is  practically  general. 

3) Scene  Independence:  The  insensitivity  of  data 
distributions  implies  nearly  scene-free  performance  of  the  optimum 
quantizer.  The  estimate  of  P  is  not  a  scene-free  quantity  in 
general.  However,  if  we  classify  scenes  into  different  groups 
according  to  the  nature  of  the  scene  (  e.g.  industrial,  residence, 
field  or  airplane,  etc)  or  the  image  points  occupancy  of  the 
scene,  then,  within  a  group,  the  amplitude  signal  can  be 
considered  as  stationary  and  the  scene  dependence  of  P  can  be 
eliminated.  Under  this  condition,  the  theory  can  be  consider  as 
scene-free . 

4)  Capability  of  real-time  processing:  In  many  image 
processing  situations,  real-time  processing  is  preferred.  The 
optimum  design  requires  knowledge  of  the  statistical  properties  of 
the  input  data  which,  in  principle,  requires  that  data  be  acquired 
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and  analysed.  Such  a  process  precludes  real-time  operation. 
However,  the  nearly  scene-free  performance  observed  in  this  study 
allows  the  designer  to  avoid  the  distribution-estimation  steps 
and  to  base  the  design  upon  assumptions  of  uniform  phase  and 
Rayleigh  amplitude  distributions.  By  making  these  priori 
assumptions  regarding  the  data,  the  quantizing  operation  becomes 
real-time . 


4.2  Distortion  Analysis 

Knowledge  about  distortion  so  far  is  limited  to  knowing  how 
much  percent  of  the  energy  belongs  to  distortion.  It  is  not 
sufficient.  In  order  to  determine  the  tolerance  to  distortion  in 
the  image,  one  has  to  know  how  the  error  energy  is  distributed. 
In  other  words,  is  it  randomly  and  homogeneously  distributed  over 
all  place?  or  does  it  build  up  coherently?  If  false  targets  are 
created,  how  do  they  relate  to  the  data  compression  procedure  and 
to  the  original  image,  etc.. 

Some  earlier  research  on  the  nature  of  image  distortion  from 
hardllmitlng  aperture  data  have  shown  that  the  hardllmlting 
procedure  does  produce  false  targets  (ghosts)  [4],  [13]— 1141 .  It 
is  also  found  that  the  distorted  image  is  related  to  the  original 
image  by  some  kinds  of  convolution  relationship  [4]  [15]. 

A  complete  solution  to  this  problem  is  out  of  the  scope  of 
this  work.  However,  a  procedure  is  developed  here  to  give  some 


useful  guidelines  to  approach  the  solution, 
By  assuming  that 

f  =  AeJ  =  E  A^e  Six  -  x^) 

*  is  * 


15'  i 

g  =  A ' eJ  =  E  AJe  si x  -  xi ) 


and  letting 


we  have 


KAi=  A’/Aj  ,  and  K$i= 
9 

g  =  E  A.'e  ^(x  -  x.  ) 

1 


(4.1) 

(4.2) 
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jKSi5l 

=  l  KAlAie  *<*  -  V 

j(K  j  1)3  j3i 

-  S  KAle  21  Ate  *.5(x  -  x^ 

j(K  j  1)31 

=  KAi*  <S(x  -  x£)  ] 


j$k 

CS  A.  e  k,5(x  -  x.  )  3 
k  *  * 


3(3,'-  2, ) 

=  CS  (A'  /Ai)  e  1  cS(x  -  xt)  ] 


CS  Ake  *k<S(x  -  xk)  3 


3(5'-  5)  33 

*  [  ( A'  /A)  e  3CA  e  3 


3  Hd  £ 


therefore 


where 


V  ha  *  x* 


(4.3) 


I  *  F(g),  hd=  P(Hd)  and  If»  F(f) 


V  'V'  V  IJf 


One  can  see  that  the  distortion  process  has  the  same  form  as 
a  linear  filtering  process  with  filter  transfer  function 


3(3'-  3) 


Hd  =  (A'/A)  e 


(4.4) 


except  that  the  parameters  of  Hd  are  discrete  random  variables. 
Let  the  complex  image 
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j 


I|®  2  Fl  <S(  u-u ^  ) 

where  =  sint1^)  represents  the  ith  pixel  location  in  a  given 

range  bin  and  represent  the  complex  pixel  value  at  that 

position  .  Then  the  complex  distorted  image  1^  can  be  expressed 
as 


V  v  hd 


*  t  2  F1  ^(u-u^  ]  *  hd 


which  shows  that  the  distortion  operation  is  equivalent  to  the 
procedure  that  first  shifts  h^  to  every  pixel  location  multiplied 
by  the  pixel  value  at  that  position,  and  then  adds  them  together 
over  the  entire  image  field. 

From  here,  one  can  immediately  observe  that 

1)  the  distortion  energy  will  be  spread  all  over  the  image 
plane; 

2)  the  nature  of  the  ghosts  depends  very  much  on  the  tail 
structure  of  h^. 

Therefore,  as  a  guideline  of  distribution  analysis,  the  shade 
of  h^  will  be  investigated  for  various  case  of  data  compression. 

4.2.1  Effect  of  Amplitude  Distortion 

For  amplitude  distortion  only,  we  have 

Hd  *  A’/A  (4.5) 

For  the  case  of  no  distortion,  Hd=  1.  The  "filter"  is  indeed 
an  all-pass  filter. 

As  the  number  of  quantized  level  decrease,  the  fluctuation 
range  of  ,ln  which  99  per  cent  of  the  data  points  of  are 
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located/  and  the  fluctuation  rate  (rate  of  changes  In  amplitude 
from  two  adjacent  elements)  increase/  so  that  more  tail  structures 
are  added  to  hd  in  the  image  domain  according  to  the  property  of 
Fourier  transform  (Figure  7).  The  probability  of  getting  false 
targets  is  therefore  Increased.  In  the  case  of  hardlimiting  (0 
bit  of  amplitude),  the  fluctuation  range  reach  its  maximum. 
Because  the  large  part  of  tail  structure  of  hd  is  in  small 
magnitude,  and  because  any  pixel  value  which  is  lower  than  certain 
threshold  (usually  -20dB)  will  be  considered  as  background  noise, 
the  major  effect  of  amplitude  distortion  is  to  increase  the  noise 
level  as  well  as  to  decrease  the  dynamic  range. 

4.2.2  Effect  of  Phase  Distortion 

For  phase  distortion  only,  we  have 

Hd»  exptS1  -  5) 

*  exp(-s5)  (4.6) 

-i  N 

hd(u)  -  N  2  expt-j  (£4^  +  v^u)) 

-i  N 

=  N  2  exp[ - j (V. (u)) ]  (4.7) 

i  1 

where  y^ is  the  Fourier  Kernel. 

Equation  (4.7)  can  be  considered  as  a  summation  of  unit 
vectors  with  angle  . 

Let  us  first  consider  two  extreme  cases: 

For  the  case  of  no  distortion.  H.  =  l,  h.  becomes  a  .5  function 

d  a 

as  expected. 

For  the  case  of  no-phase  distortion  (2'=0,  .5),i=  —4^), 

-1  N 

hd(u)=  N  2  exp(-j(4'i+  v^u))  =  0 
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Figure  7  .  Illustration  of  Analysis  of 

h  for  Amplitude  Distortion 
d 
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as  N  goes  to  Infinity.  all  vectors  are  canceled  out  Include 

broadside  because  $^s  so  as  f^s  are  uniformly  distributed  over  0 

•  •  • 

to  2tt  for  all  u.  It  follows  that  I  =h.*I.=0  and  I  =|I  1=0,  which 

g  d  f  g  g 

lead  to  an  undetermined  value  of  p.  For  N  being  finite,  this 
summation  might  not  be  zero.  However,  as  shown  in  Figure  S,  it 
does  result  in  a  meaningless  image. 

When  more  distortion  is  introduced  (the  number  of  bits 

decreases),  the  fluctuation  range  and  the  fluctuation  rate  of. the 

* 

projections  of  the  unit  vectors  Hd^on  both  the  real  axis  and 
imaginary  axis  increases  (Figure  8).  As  the  same  consequence  as 
for  amplitude  distortion,  more  tail  structure  are  added  onto  h^ . 
Therefore,  the  major  effect  due  to  phase  distortion  is  basically 
the  same  as  in  the  case  of  amplitude  distortion. 

The  above  analysis  provide  us  a  roughly  picture  about  how  the 
distorted  energy  will  distributed  and  give  guideline  for  further 
study  on  the  subject.  Although  the  distorted  image  I  is  defined 

as  the  magnitude  of  1^  and  hd  is  a  complex  quantity,  the  analysis 
on  the  structure  of  |hdl  will  still  be  very  valuable  to  the 
analysis  of  false  targets  distribution  because 


=  I  S  Fa  h(J(u-ui)  ! 

£  Z  I Fj |  !hd(u-ui) I 

which  means  that  the  result  from  the  analysis  on  I h^ I  gives  the 
information  upon  the  upper  bound  of  the  false  targets 
distribution.  In  other  words,  the  actual  result  will  never  be 
worse  than  the  result  from  the  analysis. 

4.3  Conclusions 

Summarizing  results  obtained  from  the  previous  sections,  we 
conclude: 

1)  Image  quality  is  measured  by  image  correlation  coefficient 
P  defined  by  (2.1).  It  can  be  expressed  in  the  aperture  domain  as 
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Figure  g  .  Illustration  of  Analysis  of 


h  for  Phase  Distortion 
d 
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the  product  of  the  lover  bound  correlation  coefficient  pq  and  the 
transfer  index  K. 

PQ  can  be  further  decomposed  into  two  independent  factor: 

the  amplitude  correlation  coefficient,  and  the  phase-dependent 

correlation  factor.  For  optimum  quantization,  the  amplitude 

correlation  coefficient  P^  is  equal  to  the  rras  value  of  the 

quantizer  output  level  aj  normalized  to  the  rms  value  of  the  input 

data,  with  the  input  ranges  chosen  such  that  the  optimum 

conditions  of  (2.6a)  and  (2.6b)  are  met.  For  the  case  of  the 

phase  being  uniformly  distributed,  the  phase-dependent  correlation 

factor  P*  is  the  sine  function  of  TT  /  N,  where  N,  is  the  number  of 
a  a  2 

output  levels.  For  Rayleigh  distributed  amplitude  and  uniform 
distributed  phase,  PQ  is  independent  with  the  statistical 
properties  (mean  or  variance)  of  the  image  data  and  uniquely 
related  to  the  entropies  of  the  quantized  signals. 

The  transfer  index  K  is  completely  specified  when  the 
distributions  of  the  aperture  data  and  the  quantizer  parameters 
are  given.  Given  number  of  quantized  levels,  K  is  insensitive  to 
the  changes  in  the  setting  of  output  levels.  It  has  also  shown 
that  K  is  asymptotically  independent  with  amplitude  distribution 
and  Insensitive  to  the  variations  of  the  phase  distribution. 
Therefore,  K  is  nearly  a  scene-free  parameter. 

2)  The  optimum  quantization  scheme,  which  maximizes  PQfor  a 
given  number  of  output  levels,  consists  of  two  independent 
quantizers,  an  amplitude  quantizer  and  a  phase  quantizer.  For 
each  quantizer,  the  output  level  from  each  quantizer-established- 
range  is  equal  to  the  conditional  mean  of  the  input  over  that 
range.  The  value  of  the  dividing  point  (endpoint)  between  any  two 
quantizer-established-ranges  is  equal  to  the  arithmetic  mean  of 
the  conditional  means  of  the  two  adjacent  ranges.  These  are  the 
two  conditions  for  optimum  quantization  { ( 2 . 6a ) - ( 2 . 6d ) )  .  This 
optimum  scheme  is  good  for  an  arbitrary  amplitude  distribution  and 
for  a  uniform  phase  distribution. 

3)  The  theory  of  Bit  Compression  for  microwave  imaging  is 
developed  by  presenting  an  exact  solution  of  the  analytic  equation 
of  image  correlation  coefficient  in  terms  of  the  entropy  of  the 


quantized  aperture  data  through  the  optimum  quantization 
procedure.  For  a  given  distribution  of  complex  aperture  data,  the 
theory  determines  the  minimum  number  of  bits  required  of  the 
quantized  signal  for  the  image  to  achieve  any  given  quality  or, 
conversely,  vhat  quality  of  image  can  be  obtained  for  any  given 
number  of  bits  that  the  quantized  signal  retains.  It  proves  that 
the  phase  Information  thoroughly  dominates  the  process.  It  is 
shown  that  hard  limiting  the  input  to  a  phased  array  while  leaving 
the  phase  unaltered  preserves  much  of  the  image  integrity. 
Amplitude-only  information  destroys  it.  Hard  limiting  plus  one 
bit  of  phase  information  sometimes  suffices.  Hard  limiting  plus 
two  bits  of  signal  phase  per  sample  often  results  in  acceptable 
imagery.  A  minimum  of  4  bits  of  information  per  complex  signal 
sample  guarantees  a  relatively  good  image  provided  that  3  bits  are 
used  to  represent  the  phase  and  1  bit  is  used  to  represent  the 
amplitude.  Eight  bits  are  required  to  ensure  a  high  fidelity 
image . 

4)  The  theory  provides  a  rule  of  thumb  for  optimum  bit 
allocation  between  amplitude  quantizer  and  phase  quantizer.  It 
Indicates  that  when  the  number  of  bits  per  data  sample  is  no  more 
than  8,  the  number  of  bits  assigned  to  phase  should  be  two  larger 
than  the  number  of  bits  assigned  to  amplitude.  For  finer 
quantization,  equal  numbers  of  bits  should  be  assigned  to 
amplitude  and  phase. 

5)  The  theory  Imposes  substantially  no  constraints  upon  the 
distributions  of  amplitude  and  phase  in  the  microwave  data.  It  is 
also  shown  that  the  optimum  quantization  design  derived  from  the 
theory  is  nearly  scene-independent  and  may  achieve  real-time 
performance . 

6) Flnally,  the  theory  provides  guidelines  to  explore  the 
nature  of  distortion  due  to  data  compression.  It  is  shown  that 
the  distorted  energy  is  spread  all  over  the  image  plane.  Its 
major  effect  is  to  Increase  the  noise  level  so  that  the  dynamic 
range  will  be  decreased. 
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